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ABSTRACT 

We evaluate the effect of a supersonic relative velocity between the baryons and dark matter on the thermal 
and density evolution of the first gas clouds at z < 50. Through a series of cosmological simulations, initialized 
at Zi = 100 with a range of relative streaming velocities and minihalo formation redshifts, we find that the 
typical streaming velocities will have little effect on the gas evolution. Once the collapse begins, the subsequent 
evolution of the gas will be nearly indistinguishable from the case of no streaming, and star formation will still 
proceed in the same way, with no change in the characteristic Pop III stellar masses. Reionization is expected 
to be dominated by halo masses of > 10 8 M Q , for which the effect of streaming should be negligible. 
Subject headings: cosmology: theory — early universe — galaxies: formation — stars: formation 



1. INTRODUCTION 

The formation of the first stars was a key event in the 
evolution of th e early universe (e.g. Barkana& Loeb 2001; 
i Bromm & Larsonl2004tlCiardi & Ferraral2005UGloveil2005i 
Brom met alJl2009t lLoebll2010h . After the emission of the 
Cosmic Microwave Background at z ~ 1000, the universe en- 
tered the 'Dark Ages,' the period when the distribution of mat- 
ter was very uniform and no luminous objects had yet formed. 
During this time, cold dark matter (DM) density perturbations 
grew to make the halos inside of which the first stars formed 
at z < 50. These stars are believed to have formed within 
M ~ 1O 6 M minihalos, where the infall of the baryons into 
the gravitational potential well of the DM-dominated mini- 
halo heated the gas suffi ciently to enable H2 - driven cooling 
and fragmentation (e.g. lHaiman et"aTJ 119961 ; iTegmark et alj 
1997; Yoshid a"etaLll2003l) . 

The initial growth of the density fluctuations after re- 
combination can be described using linear perturbation the- 
ory, which assumes that overdensities and velocity fields are 
small quantities. Similarly, cosmological simulations are ini- 
tialized at high z with small gas and DM peculiar veloci- 
ties, determined through a com bination of the A CDM model 
and Zeldovich approximation dZeldovichl 1 1 9701) . Recently, 
iTseliakhovich & Hiratal d2010l) added a complicating aspect 
to this picture by showing that at high redshift, there is a 
supersonic relative velocity between the baryons and DM. 
Whereas prior to recombination, photons and baryons are 
coupled such that the baryonic sound speed is ~ c/v3, af- 
ter recombination the sound speed drops to ~ 6 km s -1 . The 
root-mean square relative velocity, on the other hand, is much 
higher, 30 km s" 1 . The relative velocities are dominated by 
modes on the comoving scale of ~ 150 Mpc, the length scale 
of the sound horizon at recombination, and are coherent on 
s maller scales of a few Mpc. 

ITseliakhovich & Hiratal (120 101) examined how this effect al- 
ters the growth of DM structure, causing a small (~ 10%) 
suppression of the matter power spectrum for modes with 
wavenumber k ~ 200Mpc _1 . Using the Press-Schechter for- 
malism they have also found a decrease in M ~ 10 6 M mini- 
halos at high redshifts, z = 40. They furthermore find that the 



relative velocity effect yields a s cale-depend ent bias of the 
first halos. Extending upon this, iDalal et all ( 120101) analyti- 
cally studied the impact of the relative velocity on baryonic 
objects, finding that the collapse fraction will be slightly re- 
duced and that the large-scale clustering of M < 10 6 M Q mini- 
halos will be modulated on scales of ~ 100 Mpc. The same 
applies to any observable that traces minihalos, including 
the 21 cm absorption power spectrum. Tseliak hovich et all 
(2010) find similar results in a more detailed analysis. 

While these previous studies examined the large-scale ef- 
fects of the relative velocity, its direct influence on the delay of 
collapse and the evolution of gas falling into a single minihalo 
has yet to be considered. Simulations are necessary to under- 
stand how the relative streaming affects the non-linear regime 
and alters the processes involved in the collapse of minihalo 
gas. To this end, we perform a set of cosmological simula- 
tions which include these streaming motions. After the com- 
pl etion of this work, we became aw are of an a nalo gous paper 
bv lMaioetaT] (120 101) . Similar to lMaio et af] (120101) , we find 
a delay of gas collapse in early low-mass M ~ 10 5 — 10 6 M 
minihalos, but conclude that for typical streaming velocities 
this delay will be negligible by z ~ 1 0. Our work is comple- 



mentary to that of iMaio et al.l (120101) in that, while they are 
able to find a 1-20% overall suppression of the first objects, 
our factor of ~ 10 greater mass resolution allows us to see 
the subsequent collapse of the gas to high densities, revealing 
that even with relative streaming motions the thermal evolu- 
tion of primordial gas and subsequent Pop III star formation 
will be very similar to no-streaming cases following the initial 
collapse. 

2. NUMERICAL METHODOLOGY 

We carry out our investigation using GADGET, a widely- 
tested three di mensional smoothed particle hydrodynam- 
ics (S PH) code dSpringel et alj 120011; ISpringel & Hernquist 
2002). Simulations are performed in a periodic box with size 
of 100 h~ l kpc (comoving) and initialized at Zi = 100 with both 
DM and SPH gas particles. This is done in accordance with a 
ACDM cosmology with Q A = 0.7, Q m = 0.3, fi B = 0.04, and 
h = 0.7. We adopt erg = 0.9 for the fiducial normalization of 
the power spectrum, and also examine the case of erg = 1.4 
in which structure formation is accelerated and the first mini- 
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FIG. 1. — Top panels: Effective velocity v c ff = y^J+vJ of the gas (thin lines) and virial velocity V v j r of the simulated minihalo (thick red line). Top Left: 
'Standard collapse' case. Dashed line: v s j = 10 km s , dotted line: v s j = 3 km s , solid black line: no-streaming case. At each redshift v s was found by taking 
an average over the entire gas within the simulation box. c s refers to the average sound speed of all particles within the virial radius of the minihalo. Top Right: 
Early collapse case. Note that for the streaming cases, the redshift at which v c ff first falls below v v j r matches well with the point where the gas thermal evolution 
first follows that of V v j r . Bottom Panels: Evolution of the Jeans mass Mj with redshift, evaluated using v e ff in the role of the effective sound speed. Notation is 
the same as in the upper panels. Red line now shows virial mass M vu - of the minihalo. Green line is an exponential fit to the growth of the 'standard collapse' 
case minihalo. Gas collapse occurs quickly after Mj drops below M v j, . Note that the enhancement of v e ff due to the streaming velocity effectively increases Mj, 
causing the gas collapse to be delayed until M v j, can further grow. This alters the final gas collapse redshifts of each case (zcol — 14.4, 12.2, and 6.6 for the no 
streaming, moderate streaming, and fast streaming cases given 'standard collapse'; z co i = 23.6, 21.3, and 12.4 for 'early collapse'). 



halo collapses earlier. Each simulation box contains 128 DM 
particles and an equal number of SPH particles. The gas par- 
ticles each have a mass wsph = 8 M Q , so that the mass resolu- 
tion is, M res ~ 1 .5A^neigh»JsPH < 400 M Q , where A^neigh — 32 is 
the t ypical number of parti cles in the SPH smoothing kernel 
(e.g. iBate & Burkeri[l997h . This mass resolution allows us 
to follow the gas evolution to a maximum number density of 
"max = 10 4 cm" 3 . 

The chemistry, heating and cooling of the primordial gas 
is treated in a fashion very similar to prev ious studies (e.g. 
lBrornm&Loebll2004t lYoshida et al.1l2006h . We follow the 
abundance evolution of H, H + , H", H 2 , HJ, He, He + , He" 1 ^, e", 
and the deuterium species D, D + , D", HP, and HD + . We use 
the same chemical network as used in lGreif et al.l ( 120101) and 
include the same cooling terms. 

We first perform both the 'standard collapse' (eg = 0.9) and 
'early collapse' (a$ = 1.4) initializations with no streaming 
velocity added. For each of these we also perform 'moder- 
ate' and 'fast' streaming cases in which we include an initial 
streaming velocity v s j of 3 km s" 1 and 10 km s" 1 , respectively. 
The 'moderate' streamin g case represents the predict ed root 
mean square velocity (Tseliakhovi ch & Hiratal l2010b . given 
that peculiar velocities have decreased as (1 +z) since recom- 



bination and thus have declined by a factor of 10 at the point 
our simulations are initialized. Our v Si i values therefore cor- 
respond to velocities of 30 km s" 1 and 100 km s" 1 at recombi - 
nation, similar to the velocities chosen bv lMaio et al.l d2010l) . 
30 and 60 km s" 1 . 

3. RESULTS 
3.1. Delay of Gas Collapse 

The main effect of the relative streaming cases is to delay 
the collapse of the baryons into the DM halos. In the standard 
case, the collapse redshifts are z co i = 14.4, 12.2, and 6.6 for 
v S; i = 0, 3, and 10 km s" 1 (0, 30, and 100 km s" 1 at recombina- 
tion). The streaming cases correspond to delays in collapse of 
~7x 10 7 and ~5x 10 8 years, respectively. In the accelerated 
collapse case these values are z co i = 23.6, 21.3, and 12.4, cor- 
responding to delays of ~ 2 x 10 7 and ~2x 10 8 years. Thus, 
this delay is noticeable only for high initial values of v s ,i > 3 
km s" 1 , whereas at smaller values the delay is small compared 
to the Hubble time. 

We can understand the criterion for gas collapse in terms 
of the cosmological Jeans mass. In the usual no-streaming 
case, the slow infall of gas into the halos will first begin when 
the gravitational potential well of the minihalo, characterized 
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FIG. 2. — Effect of relative streaming on the minimum halo mass into 
which primordial gas can collapse. Each line represents the necessary halo 
masses for baryon collapse at a different redshift, marked in the plot. The 
diamonds represent the final halo masses found in 'standard collapse' simu- 
lations (zco\ = 14 for no streaming), and the squares represent masses from 
the 'early collapse simulations' (z co \ = 24 for no streaming). Note that the 
halo mass does not noticeably increase unless the initial streaming velocites 
are very high (> 3 km s _1 ). Also note that halos collapsing at high redshift 
are more affected by relative streaming, as the physical streaming velocities 
are higher at these early times. 

by its virial velocity V W1I , is large enough to assemble the gas, 
which occurs when V v - U > c s , where c s = JfJMa is the 

sound speed. Once this process begins, the sound speed c s 
will be coupled to V wu through adiabatic heating (see top pan- 
els of Fig. [TJ, and the density will scale with sound speed 
approximately as c\. In Fig. [1] we determined the proper- 
ties of the largest halo in our s imulation using the HOP tech- 
nique (Eisenstei n & Hutl ll998) to find the DM particle in the 
region of highest DM density. Assuming this particle marks 
the center of the halo, the extent of the halo was determined 
by finding the surrounding spherical region in which the aver- 
age DM density is 200p c , where p c is the redshift-dependent 
critical density. 

The bottom panels of Fig. Q] illustrate that the adiabatic 
phase of evolution will continue until the virial mass of the 
minihalo is greater than the Jeans mass of the gas, M vir > Mj. 
For the no-streaming case, we calculate Mj as 



the halo. To accommodate the cases with streaming, we re- 
place c s in equation (1) with v e ff. and the resulting increase of 
Mj is shown in the bottom panels of Fig.[JJ 

For any given collapse redshift, a larger M wu is therefore re- 
quired to trigger the collapse of streaming gas compared to the 
non-streaming case. In Fig. [2] we estimate for different red- 
shifts the mininum halo mass into which gas with various ini- 
tial streaming velocities can collapse. We arrive at these esti- 
mates using the following simple prescription. We first deter- 
mine for a range of redshifts the minihalo mass corresponding 
to a virial temperature of 1500 K, which serves as the mini- 
mum mass for collapse and cooling given no streaming. We fit 
a typical halo growth history using M v ; r (z) = Moe a: (the green 
line shown in Fig. [TJ, with Mo = 2 x 10 7 M Q and a ranging 
from -0.2 to —0.5. We vary a depending on the no-streaming 
case minihalo mass and the desired collapse redshift. For ev- 
ery given collapse redshift and streaming velocity, we first de- 
termine the redshift z eq where V V a(z) = v e ff( z )- We assume that 
z eq is the point where the gas switches from having properties 
of the intergalactic medium (IGM) to properties determined 
by the halo. Thus, for z > z eq the sound speed roughly fol- 
lows that of the IGM, c s jgm- Therefore, z eq can be found by 
considering 



Kir(z) = y/GM vil (z)/R,u(z) = x /c s , IG m(z) 2 + v s (z) 2 , (3) 



where 



fl vir (z)~210 



Mvi 



V 10 6 M f 



1/3 



^ ) .A:) PC 



(4) 



and f(z) is a factor of order unity with a mild dependence 
on redshift (Barkana & Loebl l2001). At z = z e q the effective 
gas velocity is thus v eq = V v i T (z eq )- After this point the thermal 
energy of the halo gas dominates over the energy of streaming 
motion, and its sound speed can be described by the halo virial 
velocity thereafter. 

Furthermore, z eq marks the last time that the gas density 
within R V1I is still that of the IGM. The density of halo gas 
when it first couples to the DM is then 



Peq = PlGMfeq) ; 



2x 10- 29 O n / ! 2 (l+z eq ) 3 gcm- 3 . 



(5) 



Note that, because z eq is lower for higher values of v s .j, p eq 
correspondingly decreases. 

Finally, as the gas infall into the halo continues for z < z eq , 
we estimate its average density to be 



M 



6 J g 3/ v/ 2 ' 



(i) 



Once the halo gains sufficient mass, and also provided that the 
H2-driven cooling time f coo i of the gas is shorter than its free- 
fall time fff, the Jeans and cooling criteria will be satisfied and 
the gas will begin the next phase of rapid collapse to higher 
densities, quickly reaching n max = 10 4 cm -3 . 

The cause for the delay in collapse of the streaming cases 
lies in the enhanced effective velocity of the gas, 



Veff : 



(2) 



where the streaming velocity decreases with redshift as 
v sOz) = v s.i/(l +z). As shown in the top panels of Fig. [JJ 
this delays the point at which the gas will begin falling into 



P(Z) ^ Peq 



Kir(z) 



v eq 



(6) 



The above equation describes how the gas density will adi- 
abatically evolve w ith thermal energy as it collapses (e.g. 
iTegmark et al.l fl997h . Recall that v eq is the effective sound 
speed when it first begins falling into the halo, and that adi- 



abatic evolution implies T oc p 



,7-1 



3 2/3 for 



an atomic gas 



with 7 = 5/3. Using c s oc T 1 ' 2 results in p oc c 3 . Finally, 
we replace c s with the virial velocity of the halo to arrive at 
the approximation in Equ. (6). The density will increase in 
this way until the gas virializes and reaches a maximum of 
200/3igm(z)- Inserting the applicable values for p and V v ; r at 
the given collapse redshift, z C ob we arrive at Mj(z co i). 

This model well reproduces the masses and collapse red- 
shifts found in the simulations (symbols in Fig. |2). As the 
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FIG. 3. — Evolution of temperature of the gas as density grows for both 'standard' and 'early' collapse cases. Black dots: No streaming case. Yellow dots: 
v s i = 10 km s case. Left: Standard case is shown at z = 14.4, while the v s j case is shown at z = 6.6. Right: Early collapse case is shown at z = 23.6 for no 
streaming, and z = 12.4 for v s j = 10 km s . There is almost no difference between the v s j = 3 km s (not shown) and no-streaming cases. 



streaming velocities increase, the gas density during initial in- 
fall decreases, thereby lowering the typical gas density in the 
halo and raising the minimum mass that will satisfy M v h- > 
Mj. For the average 3 km s" 1 streaming velocities, this mini- 
mum mass Mhaio will approximately double for z co \ = 30, but 
will almost stay the same by z co \ = 10. 

3.2. Thermal Evolution and Star Formation 

Figure [3] compares the thermal evolution for the 'standard' 
and 'early' collapse cases given no streaming with the cor- 
responding rapid streaming cases (v Sl i = 10 km s" 1 ). For the 
no streaming cases, the gas follows the canonical evolution 
of adiabatic heating as the IGM gas gradually becomes in- 
corporated into the growing minihalo. This gas heats to the 
virial temperature (~ 1000 K) of the minihalo until the H2 
fraction grows sufficiently high to allow the gas to cool to a 
minimum of ~ 200 K. As the gas temperature drops, its den- 
sity grows to approx imately 10 4 cm" 3 (see Bromm ^t alj2002t 
Yoshi daet al.l l2006). At this density the gas has reached the 
'loitering phase,' and this is the reservoir of gas from which 
Pop III stars will form. 

For the v S| j = 10 km s" 1 case, this evolution shows only mi- 
nor differences from that described above. The initial heating 
of the low-density gas occurs more quickly than the purely 
adiabatic rate, and the streaming velocity acts as a heating 
term for the low-density gas. However, once the gas gains suf- 
ficiently high temperature and H2 fraction, the gas cools and 
condenses to approximately 200 K and 10 4 cm" 3 , just as in the 
canonical case, though the minimum temperature is slightly 
lowered for the streaming cases. Subsequent star formation is 
therefore not suppressed, and should occur in the same way as 
it would in the no streaming case. Note also that Fig.[3]repre- 
sents significantly higher streaming velocities than typically 
expected. For the more representative v s j = 3 km s" 1 cases, 
the thermal evolution shows almost no difference from those 
with no streaming. This further strengthens the argument that 
relative streaming between baryons and dark matter will do 
little to modify Pop III star formation. 

Figure [4] further elucidates the effect of relative streaming 
on gas collapse and star formation. The delay of gas collapse 
until the minihalos have reached higher masses is evident in 



the bottom panel, which shows the fraction f$p of the mini- 
halo gas that is dense and star-forming, defined as the gas that 
exceeds densities of 1 cm" 3 . There is no star-forming gas in 
< 10 6 M Q halos for the highest streaming velocities. The gas 
fraction / gas , calculated as the gas mass in the halo over its 
total mass, is reduced by up to a factor of ~ 1.2 for v s = 3 km 
s" 1 , and a factor of ~ 1.8 for v Si i = 10 km s" 1 , even after the 
gas has reached high densities. However, once gas collapse 
has occured, the average H2 fraction /h 2 is very similar in all 
cases, as is the subsequent thermal evolution (Figure[3]l. 

4. DISCUSSION 

Our series of simulations show that Pop III star forma- 
tion will be essentially the same in cosmologies with rela- 
tive streaming motions between gas and DM, even in regions 
with streaming velocities much higher than average (v s > 3 
km s" 1 ). However, these regions of fast streaming will expe- 
rience a modest delay in the collapse redshift at which Pop 
III stars will first form, while in regions of typical streaming 
the delay will be min imal (> 10 7 years), in good agreement 
with iMaio et al.l d2010l) . In their work they also find similar 
reductions in halo gas fractions of up to a factor of 2, even 
given their larger box size and lower resolution. This is fur- 
thermore consistent with other recent work such as that of 
Tseliakh ovich et all (120101) . 

The effect on reionization should be similarly minimal. 
Though not yet known with certainty, recent work has sug- 
gested that the sources of reionization were dominated by 
early galaxies of virial temperatures above the hydrogen cool- 
ing threshold of 10 4 K (corresponding to masses > 10 8 M Q ), 
with a much smaller contribution from < 10 6 M W halos (e.g. 
Trac & Gnedinl2009tlTre"nti & StiaveUil2009tlMunoz & Loebl 
20101) . The relative streaming motions will do little to alter the 
infall of gas into the larger potential wells of ionizing galax- 
ies, and thus reionization should proceed virtually unaffected. 
The effect of streaming is mostly pronounced at the highest 
redshifts when the collapse fraction and the corresponding ra- 
diative effects of stars are exceedingly small. 

In conclusion, we have directly simulated the delay in col- 
lapse, and the subsequent thermal evolution of the first bary- 
onic structures under relative bulk velocities between gas and 
dark matter. Our results show that early star formation and 
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FIG. 4. — Evolution of gas properties with halo mass. Dashed lines: v s ; = 10 km s , dotted lines: v s i = 3 km s , solid lines: no streaming. Black represents 
the 'standard collapse' set of simulations, while blue represents the 'early collapse' set. Top panel: Average H2 fraction of the minihalo gas. Middle panel: Gas 
fraction f gas of the halos. Bottom panel: Fraction of minihalo gas that is star-forming (i.e., n > 1 cm -3 ). The reduced gas fraction and the delay of star formation 
for high streaming velocities are evident in the bottom panels. However, the H2 fraction converges to > 1CT 4 in each case, allowing for the thermal evolution of 
the highest density gas to be relatively unchanged even for the highest streaming velocites. 



subsequent evolution of reionization should quickly converge 
to the no-streaming case. Thus, results of previous and fu- 
ture cosmological studies concerning the formation of the first 
stars and galaxies will need only minimal modifications due 
to the relative streaming effect. 
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